Exact solution, scaling behaviour and quantum dynamics of a model of an 
atom-molecule Bose-Einstein condensate 
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We study the exact solution for a two-mode model describing coherent coupling between atomic 
and molecular Bose-Einstein condensates (BEC), in the context of the Bethe ansatz. By combining 
an asymptotic and numerical analysis, we identify the scaling behaviour of the model and determine 
the zero temperature expectation value for the coherence and average atomic occupation. The 
threshold coupling for production of the molecular BEC is identified as the point at which the 
energy gap is minimum. Our numerical results indicate a parity effect for the energy gap between 
ground and first excited state depending on whether the total atomic number is odd or even. The 
numerical calculations for the quantum dynamics reveals a smooth transition from the atomic to 
the molecular BEC. 
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After the experimental realization of a Bose-Einstein 
condensate (BEC) in dilute alkali gases jlj, many physi- 
cists started to consider the possibility of producing a 
molecular Bose-Einstein condensate from photoassocia- 
tion and/or the Feshbach resonance of an atomic Bose- 
Einstein condensate As discussed by Zoller jjj, 
this tantalizing problem is now coming toward resolu- 
tion. Donley et al. I| recently reported the creation of a 
BEC of coherent superpositions of atomic and molecular 
85 Rb states. This achievement is significant in that the 
entangled state is comprised of two chemically distinct 
components. 

In anticipation of this result, this novel area has at- 
tracted considerable attention from theoretical physicists 
H H. Drummond et al. emphasized the finite- 

dimensionality of the system and the importance of quan- 
tum fluctuation. Javanainen et al. M systematically 
analysed the efficiency of photoassociation of an atomic 
condensate into its molecular counterpart. In Refs. [0 |), 
large-amplitude coherent oscillations between an atomic 
BEC and a molecular BEC were predicted through 
the use of the Gross-Pitaevski (GP) mean-field theory 
(MFT). Others have gone beyond the GP-MFT HQ- 
Vardi et al. []l2| suggested that the large-amplitude atom- 
molecule coherent oscillations are damped by the rapid 
growth of fluctuations near the dynamically unstable 
molecular mode, which contradicts the MFT predictions. 
This has caused some disagreement regarding the atom- 
molecule conversion and the nature of coherence 0. In 
order to clarify the controversies raised by these inves- 
tigations, it is highly desirable to extract some rigorous 
and analytical results. 

The aim of this Letter is to show that a model Hamil- 
tonian to describe coherent coupling between atomic and 
molecular BEC's described by single modes is exactly 
solvable in the context of the algebraic Bethe ansatz 
p5[ . This makes it feasible to apply techniques well- 
established in the mathematical physics literature to 



study the physics behind this new phenomenon. The 
Bethe ansatz equations are analysed and solved asymp- 
totically in the limits of the stable molecular regime and 
the stable atomic regime. Numerical solutions are also 
obtained for the crossover regime. We identify a scaling 
invariance for the model and show that the exact solu- 
tion predicts a smooth transition from a quasi-periodic 
and stable regime to a coherent, large-amplitude, non- 
periodic oscillating regime. For finite particle number 
there are no stationary points in contrast to the predic- 
tion by Vardi et al. using a linearization scheme [ jl2| . 
Instead, the onset of strong entanglement as the detun- 
ing is decreased is identified as the point at which the 
energy gap to the first excited state takes the minimum 
value. 



The Hamiltonian takes the form 



H = -a) a + —(a^a^b + tfaa), 



(1) 



where and w denote the creation operators for atomic 
and molecular modes, respectively, £1 is a measure of the 
strength of the matrix element for creation and destruc- 
tion of molecules, and S is the molecular binding energy 
in the absence of coupling. A similar Hamiltonian was 
first used to describe optical second harmonic generation 
|l6f and (with additional damping terms) photon squeez- 
ing experiments in a two- mode interferometer |fl7f . It was 
recently investigated by Vardi et al. Jl2| as a model for 
atom-molecule BEC's. Note that the total atom num- 
ber operator N = h a + 2hi>, where h a — a^a, fib = Wb, 
provides a good quantum number since [H, N] = 0. 

Bethe ansatz solution. Following the procedure of the 
algebraic Bethe ansatz |l5|], we have derived the Bethe 
ansatz equations (BAE) which solve the model (|l|). For 
given quantum number M the BAE for the spectral pa- 
rameters {vi} read 
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Above, fc is a discrete variable which takes values 1/4 or 
3/4, dependent on whether the total number of particles 
is even or odd respectively. More specifically we have 
N = 2M + 2(k - 1/4). For each M there are 2M + 1 
families of solutions {v^ to For a given solution, the 
corresponding energy eigenvalue is 



E = 5M + 8(k - 1/4) - Q ^2 '' 



(3) 



Although we will not derive the BAE here, we remark 
that the model is obtained through a product of the Lax 
operators for the su(l, 1) algebra [|l8],[l9| and the Heisen- 
berg algebra p(| in the quasi-classical limit. The con- 
struction is similar to that used in the solution of the 
generalized Tavis-Cummings model from quantum optics 
JL8| , [L9| . The eigenstates are also obtained through this 
procedure. Consider the following class of states 
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,v M )=Y[C(v i )\V) 
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where C(v) = (u&t - |*) = |0) for = 1/4 and 

|\I/) = |0) for k — 3/4. In the case when the set of 
parameters {v^ satisfy the BAE then (^) are precisely 
the eigenstates of the Hamiltonian. 

Asymptotic analysis. In the limit of large \8/il\ we 
can perform an asymptotic analysis of the Bethe ansatz 
equations to determine the asymptotic form of the en- 
ergy spectrum. We choose the following ansatz for the 
Bethe roots 



Vi w 8 /CI + e, + fa fl /S 



i < to, 
i > to. 



Substituting into the BAE and solving for , fii leads to 
the following asymptotic result for the energy spectrum 
(cf. @) 

E m w 6(m + k - 1/4) 

+ n 2 (3m 2 - to + 4fcrn - 2kM ~ 2mM)/8. 

This result can be confirmed by second-order perturba- 
tion theory. The level spacings A m = E m — E m -\ are 
found to be 



20 2 (M + 2-3m-2fc)/5. 



We let S denote the ground state energy and A the 
gap to the first excited state. Employing the Hcllmann- 
Feynman theorem we can determine the asymptotic form 
of the following zero temperature correlations 
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where 9 = — ^a^atS + ti^aa*j is the coherence correlator. 
We then have for large 8/^N 1 / 2 ) > 

A 8 QN 1 / 2 (n„) 9 



flN 1 / 2 flN 1 / 2 8 ' N ' TV 3 / 2 
while for large S/iVlN 1 ' 2 ) < 

A 8 JIN 1 / 2 (n a ) _ fl 2 N 

9 nN 1 / 2 
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The above results suggest that the model has scale invari- 
ance with respect to the single variable QN 1 / 2 /S. They 
also suggest that the scaled gap A/^N 1 / 2 ) will have a 
minimum at some positive value of QN 1 / 2 /8 of the order 
of unity. However, as the above is an asymptotic result, 
we need to undertake a numerical analysis to obtain a 
more precise picture for the region of small ^/(flN 1 ' 2 )], 
and to establish that the scaling is also valid in this re- 
gion. 

Numerical analysis. There is a convenient method to 
determine the energy spectrum without solving the BAE. 
This is achieved by resorting to the functional Bethe 
ansatz ]22[ ]. Let us introduce the polynomial function 
whose zeros are the roots of the BAE; viz. 
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g(u)=h(i-u/i 



It can be shown from the BAE that G satisfies the dif- 
ferential equation 

uG" - (u 2 -8u- 2k)G' + (Mu -E + 8(k- 1/4))G = 

subject to the inital conditions G(0) = 1, G"(0) = 
E- 8{k- l/4)/2fc. By setting G(u) = J2n=o 9nU n the 
recurrence relation 



_ E-8(n + k- 1/4) 
— -pr, — — ttt — , , ; 9n 



n-M-1 



n(ra+l)(n + 2fc) " (n + l)(n + 2fc) 



9n-l 



is readily obtained. It is clear from this relation that g n 
is a polynomial in E of degree n. We also know that 
G is a polynomial function of degree M and so we must 
have gM+i = 0. The M + 1 roots of gn+i are precisely 
the energy levels E m . Moreover, the eigenstates (^) are 
expressible as (up to overall normalisation) 

M 

\v u ...,v M ) = Y,9n{H M - n \tftf /V) n \*) ■ (5) 

n=0 
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We have implemented these results in order to numer- 
ically solve the energy spectrum for various ranges of the 
coupling parameters. The numerical results show that 
good scaling behaviour holds for the entire range of the 
couplings, even for small particle numbers. A typical ex- 
ample for the coherence correlator is shown below. We 
remark that the result for the atomic occupation shown 
in the inset is in qualitative agreement with Fig. 1 of 
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FIG. 1. Scaling behaviour of the coherence correlator for 
the ground state as a function of the scaled detuning S. The 
curves shown are for a total atom number of N = 20, 30 
and 40. The threshold coupling for formation of the predom- 
inantly molecular BEC state is S/(Q^/N) ~ 1.4, indicating 
that there is a wide range of the scaled detuning S below 
threshold for which the ground state consists of a coherent 
superposition of the atomic and molecular states. The inset 
shows the average fractional occupation of the atomic state 
in the ground state. 



Some striking features which we have observed are that 
for fixed N there are no level crossings in the energy spec- 
trum over the entire range of couplings, and the existence 
of a parity effect for the size of the gap dependent on 
whether the total atomic number N is odd or even, as 
illustrated in Fig. 2. 

The analysis of Vardi et al. in the limit of large N 
has shown, for an entire population of molecular modes, 
the existence of two stationary points |<5/f2| — V2iV, and 
for \S/Q\ < \/2N this state becomes unstable. Our nu- 
merical analysis for finite N shows that a similar situa- 
tion occurs but the gap never vanishes. From Fig. 2 we 
see that the minimum value for A/ 5 occurs very close 
to S = nV2N. For 5/ft > V2N the state consisting 
entirely of molecules is approximately the ground state 
which is stable due to the large gap to the next excited 
level. (The same argument also applies to the negative 
S/fl, which simply follows from the inherent symmetry 
in the Hamiltonian (|l|)). 




FIG. 2. Reduction in the energy gap to the first excited 
state near the threshold coupling for formation of the predom- 
inantly molecular BEC. In other words, near threshold quan- 
tum fluctuations lead to a significant reduction in the binding 
energy of the molecules. However, note that the energy gap 
never vanishes for finite N, in contrast to the treatment of 
Vardi et al. Jl2| . The inset shows a significant parity effect. 
Near threshold the gap is larger for an odd number of atoms 
than for an even number of atoms. The solid (dash) curves 
shown are for N = 70, 90, 110 (71, 91, 111). 

Quantum dynamics. For an initial state |0(O)), the 
time evolution is given by \(f)(t)) = U{t) |</>(0)), with 

U(t) — X)m=o l m ) ( m l ex P( — iE rn t) where |m) is the 
eigenstate with energy E m . We also let g n [xn) denote the 
co-efficients in (||) for \m). Setting s = (26^6— a)a)/N, 
and taking the initial state as the pure molecular state; 
i.e, |<£(0)) = (l/%/M!)(t>) M |^), then we have 

(s(t)) = 1 - (4k- l)/N 



-jr ^ g (m)g Q (m') sin 2 (A mm -t/2) 



M 



^(n + k)g n (m)g n (m'), 



(6) 



n=0 



with 



2- n g n ^J(M - n)\(2n + 2k - ±)! 

ErLo 2" 2m (M - m)!(2m + 2k - 
E m — E m r . Similarly, if the initial state is 



and A« 



a pure atomic state, then 



s(t)) = -1 - — 22 gM(m)g M (m') sin 2 (A mm d/2) 



M 



}^M- n)g n (m)g n (m'). 



From the expression for (s(£)), we see that the short 
time behavior is quadratic rather than linear in i, due 
to the square of the sine functions, which is qualitatively 
consistent with the results of Vardi et al. |12| using a lin- 
earized model. As for the long time behavior, one can 
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see that (i) when N = 2 and 3, there are only two levels, 
so it is always periodic, (ii) when <5/f2 is large and pos- 
itive, the ground state and first excited state dominate 
the dynamics, since in this regime, the molecular state 
is approximately the ground state. As the energy lev- 
els are almost equidistant, the evolution is quasi-periodic 
(with period T « 2tt/A) and stable. A similar situa- 
tion prevails for largely negative £/f2. (hi) when S/Q is 
small, all the eigenstates are involved in the evolution 
due to strong entanglement. This means, all g n {Tn) are 
of the same order for all \m). In this regime, the levels 
are not uniformly distributed, so (s(t)) is a non-trivial 
sum of functions sin 2 (A mm /£) with coefficients depend- 
ing on g n (m) and g n (m'). Therefore, the evolution is not 
quasi-periodic. 
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FIG. 3. Time evolution of the average relative atom num- 
ber (s(t)), where r = Q.\fNt is the scaled time. The up- 
per curve depicts evolution of an initial molecular state for 
N — 40 and |<5/f2| = 6 (below threshold coupling). Clearly 
the dynamics are non-periodic. The lower curve shows col- 
lapse and revival behaviour for an initial atomic state for the 
same parameters. 
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FIG. 4. Time evolution of the average relative atom num- 
ber for iV=40 and \S/Q\ = 10 (above threshold coupling). 
The upper curve for an initial molecular state shows sta- 
ble quasi-periodic behaviour. The lower curve for an initial 
atomic state illustrates the transition from large amplitude 
collapse and revival to more stable, small amplitude oscilla- 
tions, as the detuning 5 is increased. 

Conclusion. We have shown that a two-mode model for 
an atom-molecule BEC is exactly solvable via the Bcthc 
ansatz. We have conducted both asymptotic and numer- 
ical analysis to establish scaling invariance of the model, 
identified a parity effect in the energy spectrum and in- 
vestigated the quantum dynamics. Our results indicate 
that the transition between the atomic and molecular 
BEC regimes is smooth, in contrast with Ref. JL2|. We 
believe that the exact solution will help to clarify many 
facets of the coherence properties of this model, in order 
to qualitatively compare with current experimental work 
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Vardi et al. [O claim that in the large N limit the 
molecular state is stationary when S/fl — >/2N. Observe 
however that the molecular state is never an eigenstate of 
the Hamiltonian for finite S/fl. Hence for the molecular 
state to be stationary it is implied that the gap between 
the ground and first excited state closes. For a large but 
finite N, the gap does not close (compare Fig. 2). Our 
conclusion is that there is a smooth transition from the 
large amplitude, non-periodic coherent oscillations (small 
<5/f2) to the quasi-periodic and stable regime (large S/tt), 
supported by Figs. 3 and 4. Our results are in contrast 
to Fig. 2 in Ref. Jl2| where the relative population was 
shown for a sufficiently short time interval that only the 
collapse was seen and on this basis it was mistakenly sug- 
gested that this model can describe decoherence effects. 
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